Depinning of a domain wall in the 2d random-field Ising model 
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We report studies of the behaviour of a single driven domain wall in the 2-dimensional non- 
equilibrium zero temperature random-field Ising model, closely above the depinning threshold. It is 
found that even for very weak disorder, the domain wall moves through the system in percolative 
fashion. At depinning, the fraction of spins that are flipped by the proceeding avalanche vanishes 
with the same exponent /3 — 5/36 as the infinite percolation cluster in percolation theory. With 
decreasing disorder strength, however, the size of the critical region decreases. Our numerical 
simulation data appear to reflect a crossover behaviour to an exponent /3' = at zero disorder 
strength. The conclusions of this paper strongly rely on analytical arguments. A scaling theory in 
terms of the disorder strength and the magnetic field is presented that gives the values of all critical 
exponent except for one, the value of which is estimated from scaling arguments. 



I. INTRODUCTION 

The random-field Lenz-Ising model is one of the sim- 
plest examples of random media, with applications far 
beyond magnetic systems. (For a recent review, see 
jjj.) Recently interesting nonequilibrium aspects have 
been studied, such as the field driven motion of a single 
domain wall domain coarsening in rapidly cooled 

magnetic systems, where the domain wall motion is cur- 
vature driven |f0[ , and hysteresis, with many interacting 
driven domain walls pl[ . Studies of a single driven do- 
main wall were found to describe fluid invasion in porous 
media These studies fall into the large class of in- 
terface depinning problems, which also include charge 
density wave depinning, contact line depinning, earth- 
quakes, and domain wall depinning in magnets [12 13 1. In 



these systems, second order dynamical phase transitions 
were found as the driving force F surpasses some critical 
threshold value F c , at which the interface becomes de- 
pinned, and starts to propagate through the system at a 
velocity v ~ (F — F c )^ with a critical exponent 4>. 

In previous work on a single field driven domain wall in 
the random field Ising model in three dimensions, three 
different modes of interface propagation close to the de- 
pinning threshold were identified : For weak, bounded 
disorder, the marginally stable interface at F c is facetted, 
for intermediate disorder it is self-affine, and for large dis- 
order it is self-similar. In the "faceted growth" |||],f7|]9) 
the interface propagates just as in the absence of disorder, 
penetrating the medium completely, with a roughness ex- 
ponent £ = at the depinning threshold. This type of 
interface motion can occur in any dimension, but only 
for a narrow, bounded distribution of random fields, not 
for unbounded (Gaussian) distributions of random fields. 
The existence of this phase is lattice dependent || . 

In the self-affine regime, which was seen in 3- 
dimensional simulations with bounded and Gaussian dis- 
tribution of random fields @Jl^l, neighbouring interface 
segments proceed coherently, and the interface has a 
roughness exponent smaller than one. Overhangs and 



"bubbles" (i.e., uninvaded domains left behind by the 
advancing interface) occur only below a certain length 
scale and can therefore be neglected on long length scales, 
where the interface can be described by a single-valued 
function. The critical properties of the interface near 
the depinning threshold were derived analytically start- 
ing from a continuum model with a single- valued function 
for the interface, and performing a renormalization group 
calculation, and e expansion around the upper critical 
dimension which is 5 0,0). This e expansion yields the 
roughness exponent £ = (5 — d)/Z for a d — I-dimensional 
domain wall in a in a d-dimensional system, which is 
argued to be exact to all orders in perturbation theory 
P|, however, numerical simulations show deviations from 
this prediction B. If the renormalization group result 
( — (5 — d)/3 is indeed exact, it implies that d = 2 is the 
lower critical dimension, where the ansatz of an interface 
without overhangs at large scales breaks down and con- 
ventional correlated depinning does not occur any more. 
Anisotropies in the medium may give rise to further de- 
pinning universality classes |l5f| . 

For strong disorder, the invading phase advances in a 
percolation- like manner, following routes of particularly 
high random field values. When the driving force is at the 
depinning threshold, the invading phase penetrates only 
a vanishing volume fraction of the invaded medium, just 
as a spanning cluster in percolation theory ^ffijt . Numeri- 
cal results for the fractal dimension of the invaded volume 
and the external hull of the interface suggest that this 
system is in the same universality class as uncorrelated 
site percolation dJ^]. When the disorder strength is de- 
creased, the percolation pattern coarsens, and the thick- 
ness of the percolation fingers increases and diverges at 
the critical disorder R c which marks a transition to con- 
ventional (coherent, self-affine) depinning |§,0J|]. While 
in three dimensions R c > 0, in 2 dimensions simulation 
results seem to indicate that this divergence of the finger 
width occurs only in the limit of zero disorder (R c = 0) 
||||) suggesting again that d=2 plays the role of a lower 
critical dimension. 
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This result however so far is only an indication. Within 
the numerical accuracy of the 2-dimensional simulations 
a nonzero R c could not definitely be ruled out, and the 
issue is still controversial. While 2 is the lower critical 
dimension for the equilibrium random-field Ising model, 
there is no obvious reason that this result should be tran- 
ferable to a nonequilibrium situation [ jl6| . There are 
other nonequilibrium problems with similar hurdles to es- 
tablishing the lower critical dimension. An example are 
hysteresis loops in the 2-dimensional random-field Ising 
model with many interacting interfaces. Numerical simu- 
lations, even of rather large systems (up to 30000 2 ), seem 
to converge towards zero critical disorder, but so far do 
not definitely rule out the possibility of a phase transition 
at nonvanishing disorder value either flllfl . 

Another open question is whether the type of corre- 
lated percolation found for domain-wall motion in the 
random-field Ising model does indeed belong to the same 
universality class as conventional uncorrelated site per- 
colation. In fact, while the fractal dimension of the ex- 
ternal perimeter of a site percolation cluster is 4/3 [fl7|| , 
in agreement with the dimension found for the magnetic 
interface Q , the fractal dimension of the hull (which is 
the perimeter that one measures stepping along occupied 
sites) of a percolation cluster is 7/4 M], which is different 
from the value 4/3 for the RFIM interface (see below). 

The main purpose of this paper is to present new sup- 
port for the conclusion that 2 is indeed the lower crit- 
ical dimension for the transition described above, with 
R c = 0, and to give analytical and numerical arguments 
that in 2 dimensions on sufficiently long length scales 
(even for very narrow Gaussian distribution of random 
fields) the interface propagates in the self-similar mode, 
characterized by site percolation critical exponents. In 
contrast to previous simulations of the pinned interface 
where the depinning point is approached from below, 
here it is approached from above with focus on the vol- 
ume fraction m filled by the invading phase as a function 
of the strength of the driving force (i.e., the magnetic 
field). Since m vanishes at percolation-like depinning, 
but is nonzero at correlated depinning, it can be used as 
an "order parameter" for the phase transition between 
the two behaviours. Scaling forms for m, and the width 
and length and fractal dimension of the interface, and 
other quantities, are conjectured and tested numerically, 
and ultimately justified by a scaling theory. Analytical 
arguments based on the properties of the domains of un- 
flipped spins left untouched by the proceeding interface 
give strong support for a phase transition at zero disor- 
der. They also show that the width of the critical region 
decreases towards zero as the strength of the disorder 
vanishes. Mappings of the model for several values of 
the parameters onto other percolation models with (sup- 
posedly) known values for the critical exponents confirm 
that the growth of to above depinning is characterized 
by the conventional percolation exponent (3 — 5/36. The 
numerical data are compatible with a critical point at 
zero random field, and reflect a crossover from the perco- 



lation critical exponent (3 = 5/36 very close to depinning 
to j3' = further away. The roughness exponent charac- 
terizing the interface and the fractal dimension of the in- 
terface arc C = 1 and df =4/3. Building on these results, 
we conjecture a scaling theory for the critical behaviour 
in d = 2 dimensions in the limit of small disorder. Except 
for one exponent, the exact values of all other exponents 
can be postulated, and for the remaining exponent, an 
approximate result is obtained from scaling arguments. 

The outline of the remainder of the paper is as follows: 
In section ||, the model is introduced. In sections III and 
IV, the numerical and analytical results are given. The 
final section contains a discussion of the results and ideas 
for further work. 



II. THE MODEL 

The random-field Ising model is defined by the follow- 
ing Hamiltonian 



H 



S l S ] - ^(hi + H)S l 



(1) 



The field H is the external field, denotes nearest- 

neighbour pairs, and the spin variable Si assumes the 
values ±1. In most discussions of the model, the proba- 
bility distribution of the random fields {hi} is a Gaussian 



p(hi 



exp [-h 2 /2R 2 ] 



2vri? 



(2) 



of width R. Random fields at different sites are taken to 
be uncorrelated. 

The subsequent discussion is limited to a 2-dimensional 
system on a square lattice at zero temperature. Initially 
all spins are taken to be pointing down (Si = — 1), ex- 
cept for one column of up-spins at the left boundary of 
the system, which, for convenience are given an infinitely 
positive random field hi. The upper and lower sides of 
the system are connected by periodic boundary condi- 
tions. For a given value H of the external magnetic field, 
successively all those spins are flipped up that have a pos- 
itive local field J J2j Sj + H + hi and at least one flipped 
"up" neighbour (i.e. those spins immediately neighbour- 
ing the successively propagating interface of the growing 
cluster of up spins). For negative or small positive values 
of the external field H, the interface of up spins can- 
not proceed far and stops after a small number of steps. 
For large positive H, on the other hand, the propagat- 
ing interfaces runs from left to right through the entire 
system, leaving behind only small domains of unflipped 
spins. For values H below some critical value H C (R), the 
order parameter m vanishes, i.e., m = 0, while it becomes 
finite for H > H C (R). If at R > the spin- flip avalanche 
proceeds percolation-like, the transition from to = to 
to > is continuous with 
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m~ A(R)(H - H c (R)f , 



(3) 



with some prefactor A that depends on the width of the 
random field distribution. If the invaded area has the 
same critical properties as a spanning cluster in uncorre- 
cted site percolation (in 2 dimensions), as suggested in 
@,§, then = /3p Crc = 5/36 ~ 0.139. If 2 is the lower 
critical dimension for the transition from percolation-like 
to correlated depinning, we expect that this percolation- 
like behaviour persists even for arbitrarily small random- 
ness R on sufficiently long length scales and at sufficiently 
small magnetic fields. On the other hand, if there is a 
phase transition to correlated depinning for some finite 
random field strength R Cl then the magnetization is ex- 
pected to display a jump from to some finite value m 
at at H = H C (R) for R < R c . In this paper, we will 
argue for R c = and the scenario of equation ([|) for 
R > and H sufficiently close to depinning H C (R) on 
sufficiently long length scales. If 

H + hi + (2n - 4) J > > H + h % + (2(n - 1) - 4) J , 

the local field at site i becomes positive (causing Si to 
flip up) when the nth nearest neighbour flips up. It is 
useful to define the probabilities 



(6-2n)J-H 



Pn = 



p{hi)dh t 



(4) 



(4-2n)J-H 



that a spin flips as soon as n of its nearest neighbours 
are flipped. Since we only allow spins connected to the 
advancing interface to flip, (i.e. isolated spins remain un- 
flipped), we absorb po into p%. For the advancement of 
the avalanche, there is no difference between a site that 
flips only when all four neighbours are flipped, and a site 
that does not flip even with four flipped neighbours, so we 
include in p^ all sites with hi < —H — 2 J. Subsequently, 
for convenience, we describe the system in terms of these 
four probabilities instead of H/J and R/ J, as done also 
by other authors flS|] . They span a three dimensional 
parameter space, since pi + p2 + P3 + p\ = 1- The plane 
spanned by H/J and R/J represents a cut through this 
space. Changing the external field H for a given distri- 
bution of random fields corresponds to moving along a 
line given implicitcly by Eq. (Q). There exists a critical 
surface in the parameter space that separates the region 
with m = from the region with m > 0, which contains 
the line H C (R). 

We begin by presenting some fundamental properties 
of the two limiting cases of weak and strong disorder. Fig- 
ure [j] illustrates the relation between the model parame- 
ters and the densities p n . For small R/J (weak disorder), 
regions of a given n are wide compared to R/ J, and one or 
two neighbouring values of n dominate the system, while 
for large R/J (strong disorder), the regions of given n 
are narrow, and the two boundary regions for n = 1 and 
n = 4 dominate. Obviously, in the limit R/J — * oo, the 



model corresponds to a site percolation system, where 
P4 = 1 — pi. Knowing the value of the site percolation 
threshold, pf ~ 0.59, we can immediately give an implicit 
expression for the critical magnetic field H c , 



Pi 



p(hi)dhi 



R/J 



oo. 
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FIG. 1. The densities p„ for (a) small R/J and (b) large 
R/J, represented as areas under the Gauss-distribution for 
the random field. In (a), the density pi is so small, and its 
area is so far to the left, that it is not shown in the figure. 

In the limit of small R/J, depinning occurs when H 
is such that p2 is close to one, and pi ~ y/p3- To un- 
derstand this, consider a stable phase boundary between 
the spin up and spin down regions, as drawn in figure |^. 
The invaded area (spin up) is indicated in black, as well 
as those spins that will flip as soon as one of their neigh- 
bours is flipped. The grey sites will only flip when three 
neighbours are flipped and sit therefore in the corners of 
the boundary. White sites will flip when two neighbours 
arc flipped. The depinning transition occurs when there 
exists no stable boundary that spans the system. Let us 
call a possible boundary "pinning path". The following 
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construction of a pinning path gives a good estimate for 
the relation between p\ and p 3 at the depinning thresh- 
old: Consider a system with all spins down. Now start 
at the bottom of the system at a bond that has no black 
right-hand neighbour and make a step upward. The path 
can proceed in the same direction as long as there are no 
black right-hand neighbours. It can turn left anywhere 
(if this does not lead to a black right-hand neighbour), 
but can turn right only at grey sites. Since the path 
is not allowed to intersect itself, and since it must ul- 
timately arrive at the top end of the system, the mean 
number of right turns must equal the mean number of 
left turns. Thus, for each black site that is avoided by a 
left turn, there must be a grey site, where a right turn 
can be made. Since the density of black sites is pi, the 
path has on an average 1/pi opportunities to turn left 
before encountering the next black site. After the left 
turn, the probability of encountering a grey site before 
encountering a black site, is 

Psurv = P3/(P3 + pi) ■ (5) 




FIG. 2. The phase boundary between up spins (large black 
area) and down spins (remaining area). The down spins ar 
colour coded depending on whether they flip as soon as one 
(black), two (white), or three neighbours (grey) are flipped. 

The path can survive when this probability, multiplied 
by the number of turning opportunities, is not smaller 
than 1, i.e., p^/ 1 Pi(pz + Pi) > T leading to 

P3 = Pi/(1 - Pi) 

at the depinning threshold. For small values of R (weak 
random fields), p3 and pi are also small, and ps ~ pi 
at the depinning threshold. In ||, the relation pz oc 
^i.75±o.05 wag obtained numerically (for a rectangular 
distribution of random fields), however, an exponent 2 
can also be reconciled with their figure 4, when the 
smaller slope in the lower part of the plot is ascribed 
to finite-size effects. In |LS|| , the relation p-$ ~ (2/3),o 2 
was derived, which agrees with the one given here apart 
from the prefactor. The authors of [[19[ obtained their 
relation from the condition that an initially straight in- 
terface along the first column of sites can invade the 



same number of sites in the second column of a 2-column 
system as in the third column of a three-column sys- 
tem. The prefactor should therefore change when more 
columns are taken into account. The argument presented 
in this paper does not consider the possibility of having 
two right (or left) turns in sequence, which, however, 
occurs by a factor pzj^pi + ps) less often than alternat- 
ing turns and makes therefore a negligible contribution 
to the above calculation in the limit p\ — > 0. The argu- 
ment also neglects possible correlations between different 
pinning paths starting at the same initial point. Taking 
these into account will probably change the prefactor. 
Irrespective of the prefactor, however, we can easily see 
that the density p^ is negligible in the limit of small dis- 
order: A short calculation gives the approximate result 
Ym\ R ^H c {R) ~ 2(2 - y/2)J ~ 1.172 J (also derived in 
pf), leading to 

p 3 ~ (i?/2.94J)cxp[-1.373J 2 /2i? 2 ] 

and 

Pi ~ (i?/7.95J)exp[-10.06J 2 /2i? 2 ] 

at depinning threshold. Thus, the ratio p^j p\ becomes 
arbitrarily small for small disorder, which means that 
clusters of four grey sites forming a 2 x 2-square occur far 
more often than isolated sites that do not flip with three 
flipped neighbours. Since both play the same role in the 
system by blocking an avalanche even when surrounded 
by it on three sides, and since their size difference is ir- 
relevant on long length scales, the neglection of p4 does 
not change any properties of the system in the limit of 
small disorder. 

To conclude this section, let us note that in addition 
to pinning paths that span the system and separate the 
invaded from the non invaded area, there exist pinning 
paths that are closed loops, delimiting unflipped domains 
within the invaded region. 

III. NUMERICAL RESULTS 

Earlier numerical studies |2|,|| used a quadratic system 
of L 2 sites in which the magnetic field H was increased 
incrementally, allowing the system to relax after each 
(adiabatic) increase. The system size dependent criti- 
cal field H C (R,L) at which the invading spin- up phase 
first reached the right boundary of the system, was seen 
to converge towards a constant value H C (R) as L was in- 
creased. Information about the critical field, the fractal 
dimension of the invading cluster, and its perimeter were 
obtained upon approaching the depinning threshold from 
below. 

The present work, in contrast, is mainly concerned 
with the behaviour of the fraction m of flipped spins 
above the depinning threshold H > H C (R), focussing on 
the question whether m goes to zero continuously or dis- 
continuously at depinning, and on the value of the criti- 
cal exponent (3. For this purpose, the external magnetic 
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field (or, equivalently, the densities p„) was set to a fixed 
value throughout a simulation run, and the interface was 
allowed to advance in a system of height L until it ei- 
ther came to a halt, or until it reached a cutoff distance 
which we chose to be ca. f0.5 L. Memory was allocated 
dynamically, and the system was updated by flipping all 
spins with positive local field along the interface. The 
following quantities were measured for various values of 
the system height L, and averaged over up to 250 re- 
alizations of the disorder: (i) the position of the most 
advanced and most retarded site, and the mean position 
of the interface, as well as the interface length at the mo- 
ment where it came to a stop (if it did so before running 
over the maximum allowed distance), (ii) The fraction 
of spins flipped by the avalanche (disregarding the first 
L/2 columns, where the interface had not yet reached 
its stationary behaviour, and the last columns that were 
only partially invaded by the interface), (iii) The size dis- 
tribution of the patches of unflipped spins left behind by 
the advancing interface. Before evaluation, these patches 
were allowed to relax, which is realistic for a magnetic 
system, but not for fluid invasion in a porous medium, 
where trapped regions cannot shrink. To study the low 
disorder regime, in our simulations we set p4 to zero and 
P3 to some small fixed value p 3 = 0.2, 0.1, or 0.05. p\ 
was chosen close to the depinning threshold p\(pj,). The 
snapshots in figure |^ show two pinned invasion patterns 
for two different values of the disorder (p 3 ), at pi slightly 
below the respective depinning threshold p\{pz). 





FIG. 3. The invaded area (black) for L — 400 in the low 
disorder regime (p4 = 0), slightly below depinning, with (a) 
pa = 0.2 and Pl = 0.43 < pi(0.2) ~ 0.4345 and (b) p 3 = 0.05 
and pi = 0.2236 < pf(0.05) ~ 0.226. L is the linear system 
size in the vertical direction. Roughly speaking, in the low 
disorder regime, smaller p3 corresponds to smaller disorder, 
and (pi — pf) corresponds to (H — H C (R)). 

Two characteristic trends can be discerned for decreas- 
ing p3 (i.e. decreasing random field strength R): The 
length of straight front segments increases, and the num- 
ber of unflipped domains in the invaded area decreases. 
The first feature was explained in the previous, and the 
second feature will be explained in the following section. 
First, however, let us give more details of the simulation 
results. 



A. Properties of the interface (C = 1 an d df = 4/3), 
and determination of the critical value Pi(ps) in the 
low disorder regime 




Figures ^ and [| show the probability density p(a; max , L) 
for the position a; max and a; m ; n of the most advanced 
and most retarded site of the interface for p 3 = 0.05 
and pi = 0.226 ~ pi(0.05), scaled by the system height 
L. At depinning one expects the scaling behaviour 
p(xma X ,L) ~ /i(£max/-^) with a universal roughness 
exponent £ and a universal scaling function f±. Analo- 
gously, p(x m i n , L) ~ /2(a; m in/ 'L^) with a universal scaling 
function J 2 ■ The curves for different values of L collapse 
nicely, indicating that the system is indeed at the depin- 
ning threshold, i.e., pi = p1{py,), and £ = 1. 
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FIG. 4. Probability density for the position of the most 
advanced site of the pinned interface, for p3 = 0.05 and 
pi = 0.226 ~ pi (0.05). The collapse of the curves indicates 
that C = 1. 
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FIG. 5. Probability density for the position of the most 
retarded site of the pinned interface, for pz = 0.05 and 
pi = 0.226 ~ px(0.05). The collapse indicates that f = 1. 

For larger values of pi, (i.e. H > H C (R)) the peak of 
the curves moves to the right with increasing L, and a 
nonvanishing fraction of all interfaces do not get stuck 
before reaching the cutoff distance 10.5L (especially for 
large values of L), since the system is above the depinning 
threshold. For smaller values of pi, i.e. below the depin- 
ning threshold, the percentage of interfaces that remain 
attached to the first column increases with increasing L. 
Other values of p3 give similar results. At depinning, one 
expects the mean thickness of the interface, defined as 
the number of columns spanned by the interface after it 
got stuck, to also scale as For the parameter values 
of figures H and [| this was verified with £ = 1. Once the 
interface has reached this mean asymptotic thickness of 
the order L, it becomes pinned with equal probability at 
any moment. This is reflected by the exponential tails of 
the scaling functions in figures || and ||. 



Figure |6| shows the mean length I of the pinned inter- 
face as function of L for different values of p±, at fixed 
P3 = 0.1. This length is the number of flipped spins that 
constitute the nearest- neighbour connected interface, i.e. 
the number of sites of the pinning path described in sec- 
tion [□[ Only interfaces with x m i n > L/2 and x max < 
10. hL were considered. The expected scaling form is 
/ ~ L d f g((p± — p\(p$))L 1 / v ) with universal scaling func- 
tion g, correlation length exponent u, and fractal dimen- 
sion df. At the critical threshold p\ ~ 0.317 ~ pj(0.1) 
the points do indeed follow a power law. For nearby val- 
ues of pi , the critical behaviour is only visible for L < £ 
where £ ~ ((pi — Pi(pa)) / ' P\)~ v is the correlation length. 
For larger L it crosses over to a different behaviour (linear 
in L in the case of p± > pi). The straight line in the figure 
is a power-law fit to the critical data set, with exponent 
df = 1.332±0.007. Simulations for other values of p3 give 
the same exponent df = \jv = 4/3. For a percolation 
cluster, the exponent df is 7/4 p7| . The exponent 4/3 
is retrieved, however, for a percolation cluster hull, when 
one allows steps to the next-nearest neighbours as well, 
thereby bridging the most narrow throats jL7| . Since in 
our model sites with more flipped neighbours are more 
likely to flip also, narrow throats will be bridged with 
a certain probability, thus producing the exponent 4/3. 
The same exponent df =4/3 characterizes also the frac- 
tal dimension of a self-avoiding random walk. This anal- 
ogy becomes apparent in the limit of weak randomness 
(small P3), where interfaces can be constructed by finding 
pinning paths that connect "grey" sites, as discussed in 
section [n] above. On sufficiently large scales, such paths 
are essentially self- avoiding random walks. 

100000 f — — — H 



10000 r 



1000 r 



100 




10 100 1000 10000 

L 

FIG. 6. The length of the interface / as function of L, 
for p 3 = 0.1 and Pl =0.315 (circle), 0.317 (square), 0.318 
(diamond), 0.319 (triangle up), 0.32 (plus), 0.322 (triangle 
down), and 0.325 (x). The straight line is a power law fit to 
the pi = 0.317 ~ px(O.l) data, with an exponent 1.332. 

Figure |t| shows the collapsed data of Figure [| On the 
horizontal axis, the scaling variable x = (pi — pJ(p3))L 4 / 3 
is plotted, and on the vertical axis I/L 4 / 3 , which is ex- 
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pected to be identical to the scaling function g(x) defined 
above. One can see that the scaling function is constant 
for small x, and decays with x -1 / 4 for large x. This de- 
cay corresponds to a linear dependence of I on L. The 
data points to the largest value of p\ are already outside 
the scaling regime. The scattering of the other points is 
due to not too good statistics. 
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FIG. 7. The scaling function g(x) for the interface length. 
The symbols are the same as in the previous figure. The 
straight line is a power law with exponent —1/4. 



B. The exponent (3 

Figure ^ shows the order parameter m as function of 
(pi - pi(p3))/pi for p 3 = 0.2, 0.1, and 0.05. Only those 
data points are shown that are not affected by finite-size 
effects. (When finite-size effects are present, the frac- 
tion of flipped spins decreases with increasing system 
size. The reason is that only the area behind the in- 
terface was evaluated (not taking into account the first 
L/2 columns), which cannot contain unflipped regions 
larger than L and has therefore less unflipped spins for 
smaller system size.) With the maximum system height 
L = 3200 used in the simulations, the depinning thresh- 
old could not be approached closer than shown in figure^. 
As the log-log plot shows, the data points can be fitted 
within the error bars by a power law 



m oc ((pi - pl(p?,)) / pi) 13 



(6) 



with an exponent 0.076 < (3 < 0.124 that takes decreas- 
ing values for decreasing p% and is smaller than the per- 
colation value f3perc — 5/36 ~ 0.139. Note that these 
results have to be treated with caution: (i) The data 
points cover less than a decade, (ii) Although the fitted 
power laws lie within the error bars, there appears to be 
a slight increase in slope with decreasing (pi — Pi{pz))l pi 
for all three data sets, (iii) The data are taken relatively 
far away from the critical point m = 0, in fact probably 
already outside of the scaling regime. (Attempts to ob- 
tain scaling collapses of the data using the general scaling 



form given above did not work very well.) The exponent 
/3, is actually only defined near m — as as 

[dlnm/dln((pi - pt(p 3 )) / Pi)] 

(iv) The simulated system is rather small. It is known 
from other nonequilibrium systems with quenched disor- 
der 1 11 1 that finite size effects tend to be rather large and 
result in somewhat shifted values for the critical expo- 
nents. Curiously, in |ll| collapses for the magnetization 
curves seemed also to be the hardest to obtain, (v) With 
increasing (pi — pi), the exponent (3 may cross over to 
some other value f3' . In fact, our scaling theory (see sec- 
tion IV Cj ) suggests that in the low disorder regime, very 
close to depinning, ( (p ± - pl{p?,))l Pi < {PiY^, where 



x/z is a universal exponent estimated in section [V C), 
one has m ~ {{pi - p\{pz)) I Pi) Pp "° with [3 pcic = 5/36; 
and further away ( {p\) x ^ z <C (pi — p\{pz))/ Pi )> one has 
m ^ (( (0l _ Pi{pz))l ' PiY with (3' = 0. Such a scenario 
would explain the apparent decrease of (3 with decreasing 
disorder. 
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FIG. 8. The fraction of spins flipped by the avalanche as 
function of the distance from the critical density p1{pz), for 
pa = 0.2 (solid), 0.1 (dotted), and 0.05 (dashed) for a system 
height up to L — 3200. The measured data points sit in the 
center of the error bars, and the error bars have a length of 
two times the standard deviation. The lines are a power-law 
fits to the data, with the exponents /3 = 0.124, 0.096, and 
0.076. (Note that for percolation /3 per c = 0.139.) 



Below, in section [V , we will give analytical arguments 
that for small enough p\ — p\ the exponent (3 is iden- 
tical to the percolation exponent. Furthermore, using a 
test simulation of another system for which [3 is known to 
be equal to /3 pe rc, we will show that corrections to scaling 
tend to bias the numerically observed value for (3 towards 
a "wrong" value (3 < /3 pC rc- 

Nevertheless, a few valid conclusions can be derived 
from the numerical data. First, at least for m > 0.5 
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there seems to be no tendency to approach a finite satu- 
ration value of to. If this tendency continues for smaller 
values of to, it indicates a continuous depinning transi- 
tion. Second, the critical interval p\ — p\{p^), for which 
to is smaller than some threshold value (e.g., 0.5), be- 
comes smaller with decreasing p%. Our conclusion will 
be that for sufficiently small pi — pi the exponent (3 does 
not decrease with ^3, but that the amplitude of the power 
law has to increase with decreasing p%. Below, the size 
of the critical region will be estimated using two differ- 
ent arguments, leading to a power-law divergence of the 
amplitude as P3 — ► 0. Third, using finite-size scaling, 
the fractal dimension of the invaded region can be es- 
timated. The result is compatible with the percolation 
value Df = 91/48 ~ 1.896, as found earlier in 0f 



C. The unflipped regions left behind 

Figure ^| shows the size distribution n(s) of unflipped 
domains behind the interface for ps = 0.2 and different 
values of p\ close to p\ . The size is defined as the number 
s of unflipped spins within the domain. Unflipped spins 
that are nearest or next-nearest neighbours are defined 
to belong to the same domain. (This definition allows 
for system spanning clusters of unflipped spins coexist- 
ing with system spanning clusters of flipped spins at the 
depinning point). 



l/s T f n (((pi — Pi) j ' piY n s) with f n an appropriate scal- 
ing function. This ansatz seems to give a collapse for 
r = 1.85, a n = 1.8, and pi — 0.437 (see inset). This ap- 
parent scaling behaviour, however, cannot hold for very 
small (pi — Pi)/px, since the integral sn(s)ds must be 
normalized to 1, allowing either for a scaling form with 
t = 2, or for r < 2 with no universal scaling behaviour. 
We can rule out the possibility r > 2, since the simula- 
tions as well as the analytical arguments below suggest 
that the total area of small unflipped regions decreases at 
the expense of large unflipped regions as p\ approaches 
p\. A value r > 2, in contrast, would imply that the area 
fraction covered by large unflipped regions is negligible. 

Indeed, with decreasing distance from the critical value 
PiiPi), the curve becomes flatter and does not converge 
to an asymptotic curve with finite cutoff. This behaviour 
can be interpreted as another indication that the or- 
der parameter vanishes when the critical point is ap- 
proached. If it did not vanish, the size distribution of 
unflipped domains would approach some limit distribu- 
tion with a finite cutoff at the critical point. Why this 
does not happen for a vanishing order parameter, is best 
illustrated for the site percolation limit of large disorder, 
pi = p, pa = 1 — p. There, the invaded area is identical 
(except for the first few columns) to the infinite percola- 
tion cluster. Clearly, as p is decreased towards its critical 
value p c ~ 0.59, larger and larger branches of the infi- 
nite percolation cluster become disconnected from it and 
are therefore no more flipped. All the unflipped domains 
formerly contained within this branch fuse, and larger 
unflipped domains are formed at the expense of smaller 
domains. 
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FIG. 9. Size distribution of unflipped regions for pi = 0.2 
and pi = 0.4385, 0.44, 0.442, 0.445, and 0.45 (from right to 
left). p\{p 3 = 0.2) = 0.4345 ± 0.005. The inset shows a col- 
lapse of these data, plotting s T n(s) versus ((pi — Pi)/pi) CT "s 
with r = 1.85, cr„ = 1.8, and pi = 0.437. 

Near depinning we tested the scaling ansatz n(s) ~ 



IV. ANALYTICAL RESULTS 
A. The exponent /3 

The simulation results shown in figure || give a value 
of (3, somewhere between 0.076 and 0.124, depending on 
p 3 . The obtained range does not include the percolation 
value , which is /3 pcl - c = 0.139. As argued in subsection 
IIIB, the data however are not conclusive, and a uni- 
versal value (3 — /3 pcr c (or a close value) cannot be ruled 
out. A value (3 different from /3 pC rc, would mean that any 
deviation of the densities p n from the percolation values 
should be a relevant perturbation of percolation theory. 
In other words, the conditional flipping of spins depend- 
ing on the number of flipped neighbours (i.e. P2 =/= 0, 
P3 7^ 0), should change the universality class. The case 
that (3 would depend on p% is highly unlikely. It would 
imply that even the extent to which spins are flipped as 
function of the state of their neighbours, would affect the 
value of the critical exponent. 

In the following we give several points on the critical 
surface apart from p\ = 0.59.., p± = 1 — p± where we can 
show by an analytic mapping onto percolation models 
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that (3 = /3 P erc- These results will provide a strong case 
that j3 is indeed universal. One such point is given by 

Pi = Pb, 92 =Pb(l-Pb), P3 = Pb{l-Pb) 2 ,P4 = (1-Pb) 3 , 

with pb — pi — 1/2 being the critical threshold for bond 
percolation on a square lattice. In order to understand 
this, consider a bond percolation problem, where a pair 
of neighbouring sites is connected by a bond with prob- 
ability pi,. We start with a row of up-spins at one end 
of the system and allow the avalanche to proceed to any 
site that can be accessed via open (i.e. existing) bonds. 
Some sites are invaded at the first nearest-neighbour con- 
tact with the spin-flip avalanche, others at the second 
contact, still others at the third contact, and the rest 
not even at the third contact. Obviously, the proceeding 
avalanche cannot distinguish whether it moves through a 
bond percolation system or a system with sites of differ- 
ent "colours" that are assigned according to the probabil- 
ities p n (Eq. (^)). Clearly, the invaded bonds will form a 
bond percolation cluster. Since bond percolation can be 
mapped onto site percolation on a different lattice Jl7| , 
and since the site percolation critical exponents do not 
depend on the lattice type, the number of invaded bonds 
diverges as rib ~ C(pb — pJ;) /3porc . In pc|] , it is proven 
that the number of invaded sites in a bond percolation 
problem increases with the same exponent as the num- 
ber of invaded bonds. Consequently, the order parameter 
exponent is identical to (3 perc for the above choice of the 

Pn- 
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FIG. 10. Fraction of invaded sites for bond per- 
colation with L — 3200, from a simulation with 

Pi = Pb, P2 = PbO- ~ Pb),P3 = Pb(l~ Pb) 2 , P4, = (1 - Pbf, 

with Pb ~ pi = 1/2. The solid line is a power law fit with the 
exponent (3 = 0.1296. 

Figure [Io| shows the number of invaded sites as func- 
tion of pb— pi- The simulations were performed by setting 
Pi =Pb, P2 =Pb(l~Pb),P3 =Pb(l-Pb) 2 , Pa = (1-Pb) 3 , 
with pb — pi — 1/2 in the previous simulations. The 
solid line is a power-law fit to the data, with an expo- 
nent (3 = 0.1296, which is smaller than the true exponent 
/3 P crc = 5/36 ~ 0.139 which we obtained from analyti- 



cal arguments. This suggests that corrections to scaling 
modify the asymptotic power law to an apparently dif- 
ferent power law further away from the critical point. 
Not shown in the figures is the result for the fractal di- 
mension of the interface, which is again 4/3, as in the 
previous simulations. We therefore have strong reasons 
to believe that the bond percolation case discussed in this 
subsection belongs to the same universality class as the 
simulations described in the previous section. 

For other special values of the densities p n , mappings 
on various other types of percolation are possible. For 
example, site-bond percolation (i.e., bond percolation 
where only the fraction p s of all sites are accessible) in- 
terpolates between site percolation and bond percolation 
and is obtained for the densities p\ = p s Pb, Pi — PsPb{^ — 
Pb), Pz = PsPb{l - Pb) 2 , and p 4 = p s (l - p b ) 3 + (1 - p s ). 
Mappings on short-range correlated bond percolation are 
also possible, where the probability that a bond is open 
depends on the number of bonds pointing to the same site 
(in this situation, bonds must be given an orientation). 

To summarize, there is strong evidence that the crit- 
ical exponent (3 is indeed universal and identical to its 
percolation value 5/36. 



B. The size of the critical region 



As shown in figure |8| and mentioned in section III B , 
the parameter interval during which the order param- 
eter increases from zero to a given finite value becomes 
smaller with decreasing disorder strength. Together with 
the results of the previous section, this suggests the form 
(for ( Pl - pl)/ Pl « ) 



m^Cp^(( Pl -p'i(p 3 ))/p 1 f 



(T) 



with (3 = Pp Cr c = 5/36 ~ 0.139, and with some exponent 
K. In the following, we will first derive an estimate of the 
size of the critical region and of k based on the under- 
standing that the invaded area is a coarsened percolation 
cluster, and then a different estimate that is based on a 
study of the unflipped domains. Both estimates agree 
within 12 percent, suggesting that the true value is of 
the same order as the estimated values. 



1. Arguments from percolation theory 

From the previous section, we know that the invaded 
area can be viewed as an infinite percolation cluster with 
additional sites added to it. Since a spin-flip avalanche 
can only be stopped by "grey" sites that occur with a 
probability p^ 1 , we assume now that to each site of the 
infinite percolation cluster all sites within a distance Pg 1 
are also flipped. This assumption is supported by the 
result in || that the "finger width" diverges roughly as 
p^ 1 . Clearly, when this distance becomes larger than the 
percolation correlation length £, practically all spins are 
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flipped, and the system is no more in the critical region. 
We therefore find 

at the boundary of the critical region, with v = 4/3 
known from percolation theory, leading to a size of the 
critical region 

((Pi-Pi(p 3 ))/Pi)(x^ /4 . (8) 

At the boundary of the critical region, to has some fi- 
nite value, and setting to constant in Eq. (Q), we find 

K /(3 = 3/4, or 

k= 3/3/4 = 5/48 ~ 0.104. (9) 



2. Arguments using unflipped domains 

Let us now consider the unflipped domains left behind 
by the infinite avalanche, and let us measure all distances 
in units of the "step size" / ~ l/pi • The probability that 
a (pinning) path of a given step number k has its final 
point within unit distance / from the initial point is in- 
dependent of the step size I for large I. Since there are of 
the order I 2 sites within distance I from the initial point, 
the probability that a path of k steps forms a closed loop 
vanishes as l~ 2 , which is proportional to p%. On the other 
hand, there is a cutoff 1/(1 — p SU rv/pi) to the number of 
steps of a pinning path (see Eq. (pi)), leading to 

W oc [(pi - pi(p 3 )) I 'p\Y x 

for small p 3 . As long as there exists only small and rare 
unflipped domains, the system is beyond the critical re- 
gion. The critical region can be characterized by the 
condition that the distance covered by a pinning path, 

fcmax, becomes of the same order as the distance between 
unflipped domains. Then, the picture of rare indepen- 
dent unflipped domains is no longer valid, since these 
domains can be connected by pinning paths, leading to 
a divergence of the size of unflipped domains. The above 
condition reads 

'"nix CC I CC p 3 ^ , 

and gives the following estimate of k by again putting 
to = const and using the assumption of Eq. (0) 

K= 2/3/3 = 5/54 - 0.093. (10) 

This argument also shows that the size distribution 
of unflipped domains does not converge towards a fixed 
function with finite cutoff for p\ — ► p1(pa), but rather 
that small unflipped domains become connected to form 
large unflipped domains, as observed in the simulations. 



C. Scaling theory 

The results for the order parameter and the unflipped 
domains lead to a scaling theory that is similar in spirit 
to the one proposed in |21| for the equilibrium random- 
field Ising model in two dimensions at low disorder. We 
introduce the scaling variable 

h = (pi - Pi)/ Pi, 

which measures the distance from the depinning thresh- 
old and is roughly equivalent to [H/R — H c /R), and 

t = pl, 

which is a measure for the disorder strength and is to 
leading order 

p\ ~ y/p^ oc exp[-0.343 J 2 /R 2 } . 

The characteristic length scale t^ 1 is the length of 
straight interface segments. The correlation length £ is 
most naturally identified with the diameter of the largest 
unflipped domain, weighted by the density of these do- 
mains, and can be expected to scale with some power 
of t. In equilibrium, the correlation length scales also 
as exp[CJ 2 /i? 2 ], with a constant C different from the 
depinning problem. 

Under coarse graining, the scaling variables change to 

h' = b x h t' = b z t. 

We assume that the exponent x is x = d — d/2 = 1, be- 
cause the renormalized external field H and random fields 
{hi} are given by the sum of the corresponding fields in 
the cell, leading to dimensions d and d/2 respectively 
pl[ . Here, the main assumption is that we can apply 
equilibrium rescaling under coarse graining to this essen- 
tially non-equilibrium problem. The reasoning is that in 
the limit of low disorder almost all (except for a vanish- 
ing fraction as R — > 0) coarse grained "block spins" flip 
coherently thereby optimally lowering their local energy, 
and obeying essentially the same rules as single spin flips 
on shorter length scales. The exponent z will be related 
below to the exponent k. The correlation length £ and 
the order parameter m transform under coarse graining 
according to 

£/ = £(*', h') = b~ 1 at,h) 

and 

m' = m(t',ti) = b y m(t, h) 

with y = (since all spins within a box of size b 2 are 
parallel at t = 0), leading to 

£ ~ t- 1/z l[h/t x / z ) 
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and 



m ~ t v ' z m{h/t x/z ). 



As is seen from figure [fl], on long length scales the system 
flows to the percolation fixed point at infinite disorder. 
The zero disorder fixed point and the percolation fixed 
point are connected by the depinning critical line which 
is described by the correlation critical fixed points on 
sufficiently long length scales. In the following we use 
information about the percolation fixed point to extract 
the asymptotic behaviour of the scaling functions as well. 
For h <C t x / z , we expect 

£~ (h/t x/z )- v 

with v — z^p Crc = 4/3, or 

For smaller disorder, £ is also smaller (for the same value 
of h), in agreement with our previous finding that the 
width of the critical region becomes smaller. 

From the previous section we know that for h -C t x ' z 

m{h/t x / z ) ~ {h/t x ' z f 

with j3 — Ppcrc- This gives for h -C t x l z 

m „ t to-f>*)/* h P 
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FIG. 11. Sketch of the flowdiagram for the order param- 
eter. There are 4 fixed points in this diagram: at R = 0, 
H = H C (Q) (with H c (0) = 1.172J for the model discussed 
here), which is the zero disorder fixed point discussed in this 
paper; at finite H, R — * oo, which is the percolation fixed 
point; and at H = ±oo, R — (A and B), which attract all 
the flow above respectively below the critical line. These to 
fixed points correspond to a completely flipped system and 
to a system that is not invaded at all. The thick line marked 
H C (R) is the depinning line. The arrows indicate the direc- 
tion of flow under coarse graining. The diagram explains why 
systems very close to H C (R) are dominated by percolation 
critical exponents on long length scales, as discussed in the 
paper. 



and the scaling relation 

k — {j3x — y)/2z. 

Inserting the known values for the exponents y, x, and 
/?, the relation between k and z becomes 

z = 5/72k. 

The two above estimates for k give then z = 4/7 or 
z = 9/14. 

In the opposite limit h t x / z , the order parameter 
saturates at 1, i.e., 

rh{h/t x / z ) - (h/t x / z f 

with (3' = 0. For intermediate values of h/t x ^ z , the scal- 
ing function m interpolates between the two limits. The 
exponent (3 1 is also observed for faceted growth, where 
never more than two "colours" are present, and where 
depinning occurs at p2 = 1 . In the flow diagram, faceted 
depinning occurs at the left-hand fixed point, and the 
flow follows the perpendicular axis. 

The correlation length £ is close to zero for h ^> t x / z , 
since essentially no unffipped domains are left behind. 
This means that the scaling function £ becomes propor- 
tional to {h/t x ' z )- v ' , with v 1 = oo. 



V. CONCLUSIONS 

The results of this paper confirm that 2 is the lower 
critical dimension for the transition from percolation- 
like to conventional depinning of a domain wall in the 
random-field fsing model. This conclusion is based on 
the numercial observation that the order parameter (the 
fraction of flipped spins) vanishes continuously as the de- 
pinning threshold is approached, even in the limit of very 
small disorder. Numerical results as well as analytical 
arguments show that the size distribution of unflipped 
domains left behind by the spin-flip avalanche becomes 
flatter closer to the depinning threshold, not allowing for 
a nonvanishing order parameter at the threshold. 

Furthermore, this paper supports the hypothesis that 
the order parameter exponent (5 is the same as in uncorre- 
cted site percolation. Since numerical data are not con- 
clusive and rather indicate an exponent (3 that depends 
on the disorder strength, an explicit mapping of the in- 
finite avalanche of the depinning problem onto an infi- 
nite cluster in percolation theory is performed for several 
distinct parameter values. The reason why the asymp- 
totic value of the critical exponent j3 cannot be seen in 
the simulations is that the width of the critical region 
shrinks to zero as the disorder strength vanishes. An- 
alytical arguments in this paper estimate the value of 



ff 



the exponent that characterizes the width of the criti- 
cal region, and lead to a scaling theory that relates this 
exponent to other exponents. 

As argued in ||, depinning in the random- field and 
random-bond Ising models belong to the same univer- 
sality class, and the non-equilibrium random bond and 
random field Ising model have the same symmetries near 
the critical point (for the same reasons as given in p2[), 
we therefore expect that the results of this paper are also 
valid for the random-bond Ising model. In contrast, the 
equilibrium critical behaviour of the two models is differ- 
ent. 

The question of the universality of the exponent j3 oc- 
curs also in the context of bootstrap and diffusion perco- 
lation |^3| , where all unoccupied sites of a site percolation 
problem that have a certain number of occupied neigh- 
bours are also occupied. Recently, evidence was found 
that the exponent (3 is universal in two dimensions [ p4[ . 

The scaling theory presented in this paper discusses 
only the order parameter and correlation length in the 
bulk, after the spin-flip avalanche has transversed the 
system. The scaling behaviour of the front is somewhat 
simpler for small disorder, since for small t and h the 
correlation length is identical to the length of a pinning 
path, which is proportional to 

Gfc max ) 3/4 - (th)- 3 /\ 

leading to a scaling variable th. A pinning path is essen- 
tially a self-avoiding random walk. Therefore the fractal 
dimension of the front is 4/3, which is different from the 
percolation value 7/4. However, close to the site per- 
colation fixed point (i.e. for large disorder), a crossover 
between the two exponents should be observed. Curi- 
ously, this implies that the flow diagram for the correla- 
tion length of the front shows a flow from the percola- 
tion fixed point to the t = h = fixed point, which is 
the opposite direction to the flow in figure O. A scaling 
theory for the front should therefore be performed in the 
neighbourhood of the percolation fixed point, which was 
not the focus of this paper. Other front properties like 
the size distribution of avalanches below the depinning 
threshold and the velocity of the front were not studied 
in this paper either and have still to be determined. 

It is certainly possible to generalize the scaling the- 
ory of this paper for the transition from percolation-like 
to conventional depinning to the neighbourhood of two 
dimensions by performing a 2+e expansion, in a way sim- 
ilar as in pl| | for the equilibrium model. An expansion 
around the upper critical dimension is a bigger challenge. 
Since in dimensions larger than 2 the disorder strength 
is not vanishingly small at the phase transition between 
the two different modes of depinning, the neglection of 
spontaneous spin flips away from the domain wall is re- 
alistic only under certain circumstances, for example for 
fluid invasion, for magnetic samples in a gradient field, 
or in the presence of certain long range interactions [ p2[ . 
If, on the other hand, one includes these spontaneous 



spin flips, one arrives at a hysteresis model for which an 
expansion arond the upper critical dimension 6 was suc- 
cessfully performed in |22|. 
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